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Abstract 



We discuss local duality for weak non leptonic B and A& decays in 
the rrif, — ► oo limit under the hypothesis of factorization and in the 
Shifman-Voloshin limit. We show that, under these hypotheses, 
local duality holds for the first two terms in the 1/mj, expansion 
and at the order a s in the perturbative expansion. The possi- 
ble relevance of these results for the operator product expansion 
evaluation of the 6-hadron lifetimes is discussed. 
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1 Introduction: inclusive heavy hadron de- 
cays and OPE 

In the last few years there has been an increasing theoretical interest in 
inclusive heavy hadron weak decays, mainly motivated by the possibility 
to use the theoretical method of Operator Product Expansion (OPE) to 
get information on the decay widths in the infinite heavy quark mass limit: 
m — > oo m 0. By employing OPE, a time ordered product of operators, in 
the present case the weak hamiltonian responsible for the decay, is expanded 
in terms of a sum of local operators of increasing dimension, multiplied by 
inverse powers of the heavy-quark mass. In the m — > oo limit, one may 
truncate this expansion to the first few terms, that are calculable in terms 
of perturbative QCD and some non perturbative contributions that depend 
on some known matrix elements of higher dimension operators. 

This procedure can be illustrated by taking the semileptonic inclusive 
decay width as an example. For a hadron Hj, containing one heavy quark b 
(e.g. B or At) the inclusive width is given by 

V{H h ) = -^—disc{H b \T\H b ) (1) 

ZTTIh 

where 

T = iJ d 4 xT(£(x)£(0)) (2) 

and £ is the effective lagrangian responsible for the transition. In the follow- 
ing we consider the transition b — > c only. For semileptonic inclusive decays 
one has: 

Csi = ^Vacry^l - >ft)M>f(l - 75)^ = ^V^J^ (3) 

where J% and Jjf are the hadronic and leptonic currents respectively. 

As we mentioned already, the operator product expansion, as applied to 
the product of operators appearing in @, gives rise to an expansion in terms 
of decreasing powers of the 6-quark mass m&; for semileptonic decays the first 
terms of this expansion have the form [l|]: 



r OPE(Tj . v \ G 2 F m\ ^ 2 



R(x) (l-^-)+R>(x)-^ 
\ mi mi 



(4) 
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where x = m 2 /m 2 and 



R{x) 



lx + 8ar 



x 



12x 2 logx 



R'{ 



X 



3 - 8x + 24x 2 - 24x s + 5x 4 + 12x 2 logx. 



(5) 
(6) 



The leading term in (Q) corresponds to the leading bb operator in the OPE, 
whereas the next-to-leading terms (of order 0(1/ m 2 )) arise from the color 
magnetic moment operator and the 6-quark kinetic energy; they are propor- 
tional to the hadronic matrix elements 



G = Z(H b {v)\b v 



nrT \iV (~< 

' ''° -p^b v \H b {v)) 



K 



-(H b (v)\b v ^b v \H b (v)). 



(7) 



(8) 



Here b v is the effective, velocity dependent heavy quark operator of the heavy 
quark effective theory (HQET), Z is a renormalization factor equal to one 
at the scale \l — m b . In the m b — > oo limit, G and K are finite. Let 
us incidentally note that, for the leading term in (|j), corresponding to the 
parton model result (i.e. (f|) without the corrections of the order 0(l/ml)), 
also the 0(a s ) corrections are available ||, i.e. 



F s i(b -> c£u e 



G 2 F m\ . . 2 
192tt3 |Kc61 



R(m 2 c /m 2 b ) + -JLA (m 2 /m 2 ) 

7T 



(9) 



where Aq(x) contains the perturbative corrections and can be found in 0. 

The idea that OPE can be applied to heavy hadron decays is based on the 
notion of quark-hadron duality, a concept which is widely used, but remains 
nonetheless rather vague. The physical meaning of duality is intuitively clear: 
because of the large b— quark mass, the energy release to the final hadrons is 
large as compared to Aqcd] therefore one expects that QCD can be applied 
and this would justify the approximation of identifying the hadronic decay 
width with the corresponding OPE expression. 

However a closer scrutiny shows that this expectation is in general false 
and can be retained only if some further assumptions are made. The reason 
why the physical and the OPE widths cannot be identical is in the different 
structure of their singularities: in one case there are multiparticle hadronic 
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final states and, therefore, multihadron thresholds, while in the other case 
one has quark and gluon production with a different set of thresholds. 

There are two possibilities for avoiding this negative conclusion. The 
first one consists in working far away from the physical thresholds; in this 
case OPE becomes necessarily true and one can therefore safely compute the 
hadronic quantities by their simpler OPE counterparts. However this local 
parton-hadron duality does not hold in semileptonic — > X c lv^ weak decays 
(neither it occurs in other high energy processes, e.g. e + e~ — > hadrons). The 
reason for this can be easily understood by looking at the differential decay 
rate dY / dq 2 dydvq (here q = pi + p u , y = Ei/m b , p% = m^, vq = vq/rrib). 
For q 2 and y fixed, the differential decay rate, in the complex vq plane, has 
a cut, corresponding to multihadron production. For OPE to be valid, one 
should be away from the cut, which is not the case. The second possibility, 
that can be applied in semileptonic H b decays, consists in considering, instead 
of the hadronic and OPE expressions, some averages, i.e. their smearing over 
a region of the order Aq C d @- The use of smeared quantities is generally 
referred to as global duality and has been adopted in || to prove the validity 
of OPE for semileptonic B and A& decays in a particular kinematical region, 
the Shifman-Voloshin (SV) limit [[J (see below). Considering the integrated 
semileptonic width, the authors in |4| show that OPE is valid as well, even 
without explicit smearing. The matching between the hadronic and the OPE 
expressions for the widths is proved to two orders in the 1/m expansion and 
to first order in a s . As already suggested in p, the validity of OPE is due 
to the fact that the integration contour can be deformed in such a way that 
the part which remains near the cut is of order Aqcd/ 771 an d therefore is 
vanishingly small in the m — > oo limit. 

The validity of OPE for non-leptonic decays is less clear. For non-leptonic 
transitions, neglecting penguin operators and four quarks bees operators, we 
have the effective lagrangian: 

c nl = ^v cb v: d [ Cl o 1 + c 2 o 2 ]. (10) 

Here c\ and c 2 are Wilson coefficients that in the leading-log approximation 
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are given by: 



C-2 




6/23 



6/23 




-12/23 N 



-12/23N 



and the operators O- are given by: 



O x = cy M (l - 75)6^(1 - 7 5 )m 
2 = C7 M (1 - 7 5 )Mrf7 M (l - 75)6 



(11) 



(12) 
(13) 



The OPE expression for the non leptonic decay width can be written in 
terms of the semileptonic T° PE as follows p]]: 



^OPE 
nl 



+ 



N c \V cb 

2cic 2 
~N7 



2 r OPE 
1 si 



c + c ) 1 + 



I 16 a> s 



+ 

7T 

2 /^2N3 



16 (1 — mrjmi) 



-G 



(14) 



In this case, however, we have a priori no reason to expect that the contribu- 
tion near the physical cut is small, because there is no external momentum 
q and no integration contour to deform. Therefore, for non leptonic decays, 
in order to justify OPE, one has to assume local duality or, alternatively, to 
prove it under some additional hypothesis. Clearly this puts the validity of 
OPE for non leptonic heavy hadron decays on a less firm ground 0; a signal 
of this might be the experimental result for the ratio of the total lifetimes of 
A b and B. Experimentally one has r(A&) /t(B) = 0.76 ± 0.06 (see e.g. H), 
whereas, from OPE, one expects T(A b )/r(B) ~ 1 (including corrections to 
the third order in the 1/m expansion |J). As an explanation of this discrep- 
ancy, a possible failure of local duality for non leptonic heavy hadron decays 



has been suggested |L0 



The aim of this paper is to investigate the validity of OPE for non leptonic 
B and At decays. We shall prove that local duality is indeed valid and that 
OPE can be proved under the following hypotheses: 
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1) one works in the SV || limit: 

m b , m c » 5m = m b - m c » A QCD (15) 

2) the weak non-leptonic amplitude is factorized in two parts, the first one 
corresponding to the transition b — > c and the second to a transition between 
light quarks. Since factorization has been proved only in the N c — > oo limit 
(N c = number of colours), we shall implicitly assume this limit, even if 
we retain some of the 0(1/N C ) corrections in order to keep track of terms 
violating factorization and/or duality. 

We shall show in the next section how these hypotheses are used to prove 
duality while in Section 3 we shall conclude this letter by a discussion of our 
results. 



2 B and non leptonic decays and duality 

Before considering the problem of local duality for non-leptonic B and 
decays, we shall briefly review the results obtained in H for the semileptonic 
case. 

These authors work in SV limit (|15|). The OPE result for semileptonic 
decays, as expressed by (|j) (with the corrections 0(a s ) in (|[)), becomes, in 
this limit: 




96 Sm 6 
5 nib 



G5m 4 
+ 64 + ... 



(16) 



The dots indicate several higher order corrections, i.e. terms 0(5m 7 /ml) or 
0(G5m 5 /mf) or 0(K5m 5 /ml); they can be found in [|J. As observed in [|J 
a correction 0{Ga s /mb) could also be added but has not been computed yet 
(note that G = when Hb = Ab because the light degrees of freedom have zero 
angular momentum). Computing now T s i(B — > D, D*lv) and r^(Ak — > A c lv) 
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in the same kinematical regime (|T5|), []one obtains the results 

G 2 



f lF ,2 



1927T 3 



+ 64 1 - 




64 p , - 96 <5M 6 

— 5M 5 

5 5 M B 



+ . . . 



with 



and 



5M = M B - M D = Sm 



1 - 



.Kb + G 



mi 



+ 0(l/ml) 



(17) 



(18) 



T h s ? d (A b - A c ^) 



G\ 

1927T 3 



c6 



7T 



64 



<5M 5 



96 <5M 6 ' 



+ 



with 



<5M = M Ah - M Ac = 5m 



+ C(l/m; 



(19) 
(20) 



Kb and X' in (0JD and (|0]) are given by (g) with H b = B and F b = A 6 
respectively. The inclusive semileptonic width for i? decay is saturated, in 
the SV limit, by decays into D and D* and, for A&, by its semileptonic decay 
into A c . 

Eqs. (|TBD, dT7| ) and ([0J) show that, for semileptonic inclusive decays, 
duality holds in the SV limit for the first two terms in l/ra& expansion and 
at the first order in a s ; the only term for which duality has not been checked 
in |4[] is the term proportional to a s G5M 4 /Mb appearing in ([TT|) , which has 
no matching in fllBT) (as we stressed already, the calculation of this correction 
to the OPE leading term is still missing). 

Let us show how the OPE inclusive non-leptonic decay rate (|i~4"D , can be 
matched by the sum over exclusive non-leptonic decay channels. Similarly to 
the case of semileptonic decays, we shall work in the SV limit and, therefore, 
we shall take the results of |1, i.e. that duality is proved for the first two 
terms in the 1 / mj expansion and at the first order in a s for the semileptonic 

^^Under the SV kinematical conditions, the form factors can be expanded around the 
zero recoil point and corrections to this limit appear only at the order 1/M 2 . 
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decays. As a second hypothesis, we shall assume factorization of the decay 
amplitudes. Because of our hypotheses, we expect that, also for non leptonic 
decays, duality may be proved only for the first two leading terms in l/m& 
expansion and up to order a s . 

The full width for non-leptonic decays is given by 

T nl = discli f d 4 x(H b (v)\T(C nl (x)C n i(0))\H b (v)) 

ZM Hb I J 



YP^ip-Px) \{X\C nl (Q)\H b (v))\ 2 (21) 



where C n i is the effective lagrangian for non leptonic decays given in (pT0|) . 
Let us now write 

|X> = \X U X C ) (22) 

where X u is a set of light particles, and X c is a set of hadronic charmed 
states. Let us now consider the matrix element 

(X u X c \C nl (0)\H b (v)} . (23) 

Assuming factorization, i.e. inserting the vacuum in all possible ways, we 
obtain, by using the Fierz identities, 

(x u x c \c nl (o)\H b (v)) = (<* + ?Q 9z Vcb v: d (x u \d^(i - 75 H0> x 

(X c \^(l- l5 )b\H b (v)). (24) 

A comment on the factorization procedure we have adopted is now in 
order. From eqs. (|lOD ,(|l2D and (|13"D , we see that, using the Fierz theorem, 
one can write for £ n i(0) the following expressions: 

Aa(0) = ^V cb V: d [( Cl + ^)0 1 + OJ (25) 

or 

Au(0) = ^WftKoa + |r)O a + 6 2 ] . (26) 

Oi and O2 are written as products of two coloured currents; in the factoriza- 
tion approximation they give no contribution because the matrix elements 
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of a coloured current between physical hadronic states vanish. In exclusive 
decays, (|25|) and (|^) produce respectively the so-called class I and class 
II decays. We shall now show that, within our hypotheses, we can always 
choose the form (|25|) for C n \. In order to be definite let us consider Hb = B° 
decays. If \X U X C ) contains D + or D* + states, we can assume eq.(Po]) and 
apply factorization as in (p4|), with \X C ) = \D + ) or \D* + ) and neglecting 0\. 
If \X U X C ) contains D° or D*° states, we can evaluate the matrix element in 
two ways. 

1) We can write £ n z(0) according to (|25|) and neglect 0\. In this case one 
shall factorize the amplitude as follows (considering only D° for simplicity): 



c nl (o)B«) = [ Cl + ^^=v cb v: d (x-\dY(i~iMo)y< 



X + D° 



cf(l- 78 )b\H h (v)). (27) 



where X + ,X~ are multiparticle states containing only light hadrons (the 
superscript ± indicates the total charge). On the basis of the results for 
semileptonic decays contained in [4| and |J, ( p7]) represents in the SV limit 
a higher order correction (of the order of 0(SM 2 /M^)) because B° decays 
semileptonically only to D + , D* + in this limit. 

2) The second way to evaluate the contribution of these states is by 
making use of fl26|) and neglecting, in the factorization approximation, 0<i- 
However one can see immediately that, also by this choice, the matrix element 
is a higher order correction and should be neglected. As a matter of fact, since 
we work in the N c — > oo limit, where factorization can be proved to be valid, 
we should evaluate the resulting expression in this limit (for discussion on this 



point see |llf , ||12|| ). But in the N c — > oo limit these states would contribute 
to the width a term C(c 2 + Ci/N c ) 2 = 0(c 2 ) 2 = 0(a 2 ); in the approximation 
when one neglects higher order 0(ot 2 s ) perturbative contributions, this term 
should be neglected as well. 

On the basis of this discussion we conclude that, in the SV limit, and 
taking into account only the first two terms in the l/m?, expansion and at 
order a s in the perturbative series, we can limit ourselves to consider only 
class I decays, as expressed by (Ell) . 
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From ( pip and ([24]) one obtains 



r*{H b ) = ^\V cb V:f ( Cl + ^) 2 / ^3^(9, ^^(g), (28) 

where is the same hadronic tensor appearing in the semileptonic b — > c 

decays 

^1 = (27r) 3 £5 4 (^-g-^ c )x 

(tf b |&7"(l - 75)c|X c )(X c |c 7 "(l - ls)b\H b ) (29) 

while T^y is given by: 

T^(g) = (2vr) 4 5> 4 (g-g x Jx 

(0|M 7M (l-75)rf|X M )(X u |J 7! ,(l-7 5 ) M |0) . (30) 

The vector part of the tensor T^ u is theoretically known for large values 
of q 2 up to the order o? s from e + e~ into hadrons: here we need only the 
expression to the first order in a s . We shall write 

T^(q)=T^(q)+T^ CD ( q ) , (31) 

where TQ CD (q) is the QCD contribution to T^ u (q); neglecting the u and d 
quark masses, it is given by: 

T% D {q) = 2N c j- (l + ^) ( ? V - q 2 g^)0(q 2 - Qo) ■ (32) 

The factor 2 with respect to the result for R e + e - is due to the contribution 
of both vector and axial currents in (|30|) ; q 2 ~ 1 — 2 GeV 2 is the onset of the 
QCD behaviour. Tj^ u (q) takes into account the contribution of the low lying 
{m 2 R < <7q) particles: it, a±, p, p'. It can be easily checked that, in the SV limit, 
these resonances contribute to the non-leptonic rate with terms of the order 
G 2 F 5m 3 fft, where fn is a mass parameter, of order Aqcd, characteristic of 
the resonance. As we shall discuss below, these terms represent corrections 
0(Aqcd/ '5m) 2 to the leading contribution arising from Tq CD and will be 
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neglected because, in the SV limit, 5m » Aqcd (see eq. (|T5|)). We shall 
therefore put 

T^(q)^T& D (q) . (33) 

To obtain the non leptonic decay width T n i(H b ) from (^) we observe that 
for the semileptonic inclusive b —>■ c decay, the following formula holds 

r s i(H b ) = ^f\V cb \ 2 J ^W%{q,v)T%,{q) , (34) 

where, for q 2 > q 2 , the leptonic tensor T l differs [] from (^) only by the 
overall factor N c (l + a s /n): 

TlM) = ^{q»q v -q 2 g» v ) ■ (35) 

From the previous formulae one obtains an expression for the non leptonic 
width, computed in the factorization approximation, as a function of the 
semileptonic one: 

T h n f(H b ) =N e (l + ^) \V cb \ 2 ( Cl + §-J 2 T h rm + ... (36) 

where the omitted terms are O (A 2 / 5m 2 ) or 0{5m 2 /m 2 ). In the SV limit 
we can substitute in this equation expressions ([17]) and (|19D for B and A& 
respectively. 

To compare this result with (0), we develop the Wilson coefficients C\ 
and C2 to the order a s : 

O.00 = l--|^^log(^/A4-) 

cM = 3^Mlog(^ 2 /A/a-) • (37) 



We notice that the combination c\ + 02/ N c in (p6|), does not contain first 
order a s corrections. By developing ( |T4] ) in the SV limit, we obtain, at the 

2 The presence of q$ 7^ in ( |32| ) introduces corrections 0(Aq CD / '5m 2 ), as it can be 
easily verified. Similarly to our discussion above, we neglect this correction in the SV 
limit. 
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first order in a,: 



^OPE 
nl 



GpN, 



1927T 3 



'cb\ 



64 X 5 

— dm 

5 



96 5m 6 G5m A 

+ 64 + 



5 m b 



m b 



20cic 2 G 
N r 5m 2 



(38) 



On the other hand, from 
the expression 



6|), after substitution of ( pf ) or ( |I~9"D we have 



G F N C 

1927T 3 



eft 



64 p „ 96 5M 6 

— 5M 5 h 64 h . . . 

5 5 M Hh M Hh 



(39) 



which is valid in the same limit. Comparing these two equations, we see 
that, a part from the term 2OC1C2G / (N c d~m 2 ) , which is present in (|38D and 
has no counterpart in (|39D , there is matching between the two expressions. 
We notice that, in this term, G is of the order A.q CD ; therefore it represents 
a correction of the order a s AQ CD 5m 3 and therefore belongs to a class of 
corrections that we have neglected (for example those coming from the low- 
lying resonances)^. Moreover it represents a correction 0(1/N C ) and should 
be neglected, in the approximation we are considering. 

The previous analysis shows that the inclusive non leptonic B and 
decay widths, computed in the Shifman-Voloshin limit (|I5| ) by factorization, 
are equal to their OPE expressions including the leading term O (5m 5 ) and 
the next-to- leading term 0(5m 6 /m%) in the OPE, at the order a s in the 
perturbative expansion and up to terms 0(AQ CD /5m 2 ) that are negligible in 
the SV limit. Within this limitations our analysis shows that local duality is 
indeed satisfied by non leptonic heavy hadron decays. 



3 Conclusions 

Let us summarize and discuss our results. We have proved local duality for 
B and A b non leptonic decays, i.e. the equalities 

T°T E (H b ) = T h n f(A b ) = T h n f(B) , (40) 

3 In this approximation one should neglect the term 64G5m 4 /nib as well. It can be 
noted, however, that this term has an exact counterpart in the hadronic expression. 
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for the first two terms in the l/m& expansion and at the order a s in the 
perturbative expansion. These results have been proved in the Shifman- 
Voloshin limit and assuming factorization of the weak amplitudes. 

One may wonder on the experimental significance of SV limit; for exam- 
ple, in this kinematical regime, the B meson would decay semileptonically 
only to D and D*. This prediction however, does not agree with the data, 
since experimentally one has 

T(B^D,D*£v e ) , . 

r( B ^W -°- 64± °- 08 (41) 

which is obtained considering both B° and B^ semileptonic decays |0|]. 
This discrepancy is expected since the assumption fll5]) m c » 5m is clearly 
violated. 

Secondly, our results are based on the assumption of factorization of the 
hadronic final state into two parts, one containing a single charmed particle 
(D, D* or A c ) and the other containing only light hadrons. Factorization 
holds in the N c — > oo limit, and can therefore be violated by terms O(l/N ). 
Moreover, there are arguments that can justify factorization in the m& — > oo 



limit []I4 , but it is not clear to which order in the l/m& expansion it holds. 
In [15 1, an estimate is given of the nonfactorizable part of the exclusive pro- 
cess B° — > D + 7i~: they found a term proportional to the matrix element G 
(see eq. ([?])) and color-suppressed with respect to the leading one; moreover 
these authors find that in the SV limit such a term is not suppressed by 
powers of 1/mb. This result agrees with our eq. (|38|) , where similar terms 
are present (those proportional to 2OC1C2G / \N c 8m 2 )) . As discussed above, 
besides being colour suppressed, in the SV limit these contributions are neg- 
ligible because are of order Ag C . D /5m 2 ; it is quite possible that, beyond the 
SV limit, nonfactorizable terms arise that are not suppressed by 1 jm\ in 
the infinite heavy quark mass limit. In [[U| the problem of factorization in 
inclusive non leptonic B decays is discussed assuming a standpoint comple- 
mentary to the one taken in the present paper: these authors assume local 
duality from the outset and discuss the limits of factorization. As far as our 
approach is comparable with theirs, the results of the present paper agree 
with those of fl6| . 
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Finally let us note that, because of the assumed SV limit, the non leptonic 
widths depend on SM and not on Mn b . Because of (|18|) and (p0|) , both for B 
and Ah decays 5M = 5m + 0(l/ml) and there is no way for a OiXjm^) term 
to appear in the ratio r(Ab)/r(B), contrarily to the possibility suggested in 
[|IIJ. This result, however, might be heavily dependent on the kinematics of 
the SV limit and the possibility remains that, beyond this limit, 0(l/mb) 
corrections to the non leptonic B and lifetimes appear. 
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